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@ Verfahren zur Erhohung der Leistungsfahigkeit einer Computereinrichtung bei Finite-Elemente-Simulationen 
und eine solche Computereinrichtung 

@ Verfahren zur Erhohung der Leistungsfahigkeit einer 
Computereinrichtung zur Finiten Elemente-Simulation 
durch automatische Generierung geeigneter Basisfunk- 
tionen mit Hilfe von uniformen Tensorprodukt B-Splines 
mit den folgenden Schritten: 

- Definieren (1) eines Simulationsgebietes (Q) und Ab- 
speichern von Gebietsdaten des Simulationsgebietes (Q), 

- Eingabe (2) und Abspeichern von Randbedingungen, 

- Festlegen (3) einer vorgebbaren Gitterweite (h) und ei- 
nes vorgebbaren Grades (n) der B-Splines, 

- Ermitteln (4) eines sich uber das Simuiationsgebiet (H) 
erstreckenden Gitters und eines Typs der Gitterzellen, 

- Klassifizieren (5) der B-Splines, deren Trager das Simu- 
iationsgebiet (Q) schneiden, durch Ermitteln von inneren 
und aufceren B-Splines, wobei fur innere B-Splines der 
Schnitt des Tragers mit dem Simuiationsgebiet (Q) ein 
Quadrat bzw. einen Wurfel einer vorgebbaren Kantenlan- 
ge (s), welche mit der Weite des Gitters (h) ubereinstim- 
men kann, enthalt, 

- Ermitteln (6) von Kopplungskoeffizienten <Oj s) zur Bil- 
dung von Linearkombinationen von inneren una aufceren 
B-Splines durch Lagrange-lnterpolation auf zu den auRe- 
ren Indizes nachstliegenden Feldern innerer Indizes, 

- Speicherung und Ausgabe der die Basisfunktionen be- 
stimmenden Parameter. 
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Beschreibung 

[0001] Die Erfindung betrifft ein Verfahren zur ErhShung der Leistungsfahigkeit einer Computereinrichtung bei Finite- 
Elemente-Simulalionen durch effiziente automatische Konstruktion geeigneter Basisfunktionen fur die Berechnung von 

5 Naherungslosungcn und cine solche Computereinrichtung. 

[0002] Eine Vielzahl technisch-physikaUscher Phanomene lasst sich durch partielle Differenzialgleichungen beschrei- 
ben. Hierunter fallen u. a. Probleme aus der Stromungsmechanik (z. B. Stromung um eine Tragflache), der elektroma- 
gnetischen Feldtheorie (z. B. elektrischer Feldverlauf in einem Transistor) oder der Elastizitalstheorie (z. B. Deformation 
einer Autokarosserie). Die genaue Kenntnis und Beschreibung derartiger Vorgange ist ein zentraler Bestandteil bei der 

10 Konstruktion und Optimierung technischer Objekte. Um zeit- und kostenintensive Experimente einzusparen ist das In- 
teresse an computergestutzten Simulationen groB. Hierbei haben sich Finite-Elemente- Verfahren (FE- Verfahren) eta- 
bliert und sind seit langem Gegenstand intensiver Forschung. Dies betrifft insbesondere auch automatische Netzgenerie- 
rungsverfahren als Grundlage flir die Konstruktion geeigneter Basisfunktionen. 

[0003] Fig. 1 veranschaulicht in der linken Bildhalfte den Stand der Technik beim Ablauf einer FE-Simulation fur ein 
15 lincares Randwertproblem als lypischcs Modellbcispicl. Ausgehend von den Daten, die die Geometric des zu simulic- 
renden technischen Objektes beschreiben, wird zunachst ein System von Basisfunktionen konstruiert, das einerseits die 
Einhaltung der Randbedingungen ermoglicht und andererseits zur Approximation der gesuchten Losung geeignet ist. 
Dann wird mit Hilfe dieser Basisfunktionen ein lineares Gleichungssystem mittels numerischer Integrationsverfahren 
aufgestellt. SchlieBlich werden die Koeffizienten der gesuchten Approximation als Losung dieses linearen Gleichungs- 
20 systems crmittelt. 

[0004] FE- Verfahren oder deren Anwendung sind Gegenstand einer Reihe von Patenten. Beispiclhaft sei US 4,819,161 
genannt, wo die FE-Approximationen einer groBen Klasse von Differenzialgleichungen automatisiert werden, sowie 
US 5,731,817, wo ein Verfahren zur Generierung von Hexaeder-Netzen die Grundlage fur ein Verfahren zur FE-Simula- 
tion bildet. 

25 [0005] Bei den meisten praxisrelevanten FE- Verfahren werden Basisfunktionen verwendet, die auf einer durch Vernet- 
zung des Simulationsgebietes erzeugten Zcrlegung definiert sind. Fig. 2a zeigt cine Auswahl gebrauchlicher Elemente 
deren Dimension, Grad, Glattheit und Parameter in Fig. 2b aufgelistet sind. Eine Obersicht iiber Verne tzungsmethoden 
ebener Gebiete findet sich beispielsweise in K. Ho-Le, Finite element mesh generation methods: a review and classifica- 
tion, Comp. Aided Design 20 (1988), 27-38. Die Vemetzung komplizierterer dreidimensionaler Gebiete gestaltet sich 

30 nach derzeitigem Kenntnisstand uberaus schwierig, wie S. Owen, A survey of unstructured mesh generation technology, 
Proceedings, 7th International Meshing Round T\ible, Sandia National Lab (1998), 239-257, zeigt. Die Verfahren beno- 
tigen lange Rechenzeiten und sind teilweise auch noch nicht vollkommen automatisiert. Es gibt jedoch in jiingerer Zeit 
eine Reihe sehr innovativer neuer Ansatze. Beispielsweise wird in A. Fuchs, Optimierte Delaunay-TViangulierungen zur 
Vernetzung getrimmter NURBS -Korper, Universitat Stuttgart, 1999 eine Krafteverteilung simuliert, um eine optimale 

35 Verteilung von Triangulierungspunkten zu erzielen. In US 5,729,670 werden zwei- und dreidimensionale Netze durch 
Losen von Stromungsproblemen crzeugL, was eine interessante Umkehrung des ublichcn FE-Mechanismusses ist. Dane- 
ben wurden viele Algorithmen entwickelt, um einzelne Teilaspekte der Netzgenerierung zu verbessern. Beispielsweise 
wird in DE 196 21 434 Al und US 5,774,696 ein Verfahren zur Elimination von Uberschneidungen mit vorgegebenen 
Kanten oder Grenzflachen bei Delaunay-Triangulierungen beschrieben. 

40 [0006] Vernetzungsfreie FE- Verfahren haben bisher keine Bedeutung fur die Anwendungen erlangt. Sowohl bei dem 
Lagrange-Multiplier- Verfahren, s. z. B. J. H. Bramble, The Lagrange multiplier method for Dirichlet's problem, Math. 
Comp. 37 (1981), 1-11, als auch bei der Penalty- Methode, s. z. B. P. Bochev und M. Gunzburger, Finite element methods 
of least squares type, SIAM Review 40 (1998), 789-837, steilt bei Verwendung einfacher, stabiler Basisfunktionen die 
Behandlung von Randbedingungen ein Hauptproblem dar. 

45 [0007] Bei vielen technischen Simulationen ist die automatische Netzgenerierung sehr aufwendig und nimmt den weit- 
aus grbBten Teil der Rechenzeit in Anspruch. Des Weiteren ist die Approximationskraft der ublicherweise verwendeten 
linearen und rnultilinearen Basisfunktionen gering. Um genaue Resultate zu erzielen, muss also eine groBe Anzahl von 
Basisfunktionen verwendet und damit ein entsprechend groBes Gleichungssystem gelost werden. Ansatzfunktionen ho- 
heren Grades uber Triangulierungen besitzen in der Regel ebenfalls ein ungunstiges Verhaltnis zwischen der erreichba- 

50 ren Genauigkeit und der Anzahl der verwendeten Basisfunktionen. SchlieBlich lassen sich auf unstrukturierten Netzen 
keine glatten Basisfunktionen auf einfache Weise definieren. Bereits fur stetig differenzierbare Elemente sind sehr spe- 
zielle Konstruktionen notwendig (s. Fig. 2a). 

[0008] HAUSSLER-COMBE, U., u. a.: "Die elementfreie Galerkin Methode - ttbcrblick und Anwendungsbeispiele", 
Bauingenieur, Band 74, Nummer 3, Marz 1999, Seite 130 bis 138 beschreibt ein Verfahren, bei dem eine Familie von Ge- 

55 wichtsfunktionen dazu verwendet wird, Polynome zu lokalisieren und damit eine Basis zu konstruieren. Randbedingun- 
gen werden mittels des sogenannten Penalty-Ansatzes beriicksichtigt, wobei die Wahl der Konstanten erhebliche Pro- 
bleme aufwirft. Die Wahl der Gewichtsfunktion gemaB Bild 3 dieses Dokuments dient zum Modellieren eines Risses, 
und nicht - wie gemaB der vorliegenden Erfindung - zur Einhaltung von Randbedingungen. GemaB der vorlicgenden Er- 
findung hat die Gewichtsfunktion eine grundsatzlich andere Bedeutung, namlich die, dass homogene Randbedingungen 

60 per Konstruktion exakt eingehalten werden und sornit gerade der in diesem Dokument angegebene Penalty- Ansatz uber- 
fliissig wird und die damit verbundenen Probleme beseitigt werden. Die Lokalitat der Basis ist gemaB der vorliegenden 
Erfindung durch der Verwendung der B-Splines gesichert. 

[0009] B ARTELS, R. H., u. a.: "An Introduction to Splines for Use in Computer Graphics & Geometric Modeling", 
MORGAN KAUFMANN PUB,, 1987, Seite 46 bis 63, steilt einige Grundzuge der B-Spline Theorie dar. Ein konkreter 
65 Bezug zur Finite-Elemente-Methode ist aber nicht offenbart oder nahegelegt. 

[0010] V. L. Rvachev und T. I. Sheiko: "R-Functions in boundary value problems in mechanics", Applied Mechanics 
Reviews, 48 (1995), Seiten 151 bis 188 und V. Shapiro: "Theory of R-Functions and Applications: A Primer", Cornell 
Programmable Automation, Sibley School of Mechanical Engineering, Ithaca NY (1988), CPA 883, geben Verfahren an, 
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wie Gewichtsfunktionen fur Gebiete konstruiert werden konnen, die sich speziell durch Schnitt, \fereinigung und Kom- 
plementbildung aus elementaren Teilgebieten erzeugen lassen. Diese Gewichlsfunklionen konnen zur Einhaltung von 
wesentlichen Randbedingungen verwendet werden. 
[0011] V. L. Rvachev und T. I. Sheiko, V. Shapiro und I. T^ukanov: 

"On Completeness of RFM-Solution Structures", Computational Mechanics, 25 (2000), Seitcn 305 bis 3 16 und V Sha- 5 
piro und I. Tsukanov: "Meshfree Simualation of Deforming Domains", Computer-Aided Design, 31 (1999), Seiten 459 
bis 471 setzen zwar B-Splines als Basisfunktionen ein, das Stabilitatsproblem bleibt aber ungelSst. 
[0012] R Oswald: "Multilevel Solvers for elliptic problems on domains", MuUiscale Wavelet Methods for PDEs, W. 
Dahmen et al. (Herausgeber), Academic Press (1996), Seiten 3 bis 58, beschreibt ein Verfahren, wie B-Splines durch eine 
Erweiterungskonstrukuon slabilisiert werden konnen. Die dabei verwendete Konstruktion unterscheidet sich aber von 10 
der erfindungsgemaBen Konstruktion. Dariiber hinaus werden wesentliche Randbedingungen bei der Konstruktion nicht 
beriicksichtigt. 

[0013] Der Erfindung liegt das Problem zugrunde, die Leistungsfahigkeit bekannter FE- Verfahren und die FE- Verfah- 
ren ausfuhrender Computereinrichtungen durch die effiziente Konstruktion von Basisfunktionen mit gunstigen Eigen- 
schaften zu crhohen. Insbcsondere soil die Vemetzung des Simulationsgebictes vollstandig enlfallen, evcntuell gegebene 15 
Randbedingungen eingehalten werden, genaue Losungen mit relativ wenigen Koeffizienten erreichbar und das entste- 
hende Gleichungssystem effizient losbar sein, Dadurch sollen die Nachteile des Standes der Technik uberwunden und so- 
mit die Genauigkeit und Geschwindigkeit der Simulation physikalischer Eigenschaften bei der Konstruktion und Opti- 
mierung technischer Objekte verbessert werden. 

[0014] Einige zentrale BegrifTc und Bezeichnungen, die in der nachfolgenden Beschrcibung des erfindungsgemaBen 20 
Verf aniens verwendet werden, sollen vorab erlautert werden. 

[0015] Das Simulationsgebiet ist eine beschrankte Menge der Dimension d = 2 oder d = 3, auf der die zu untersu- 
chenden physikalischen GroBen mittels der FE-Methode approximiert werden sollen. Der Rand des Simulationsgebietes 
wird mit T bezeichnet. Unter einem Gitter mit Gitterweite h versteht man eine Zerlegung einer Teilmenge der Ebene oder 
des Raumes in Gitterzellen Z^. Dabei ist jede Gitterzelle je nach Dimension d ein Quadrat oder ein Wiirfel mit Kanten- 25 
lange h. Genauer gilt Zk = kh + [0, h] d , wobei k eine Menge von ganzzahligen d- Vektoren durchlauft. Die uniformen Tcn- 
sorprodukt B-Splines in d Variablen vom Grad n zur Gitterweite h werden mit bezeichnet, s. z. B. C. de Boor, A Prac- 
tical Guide to Splines, Springer, 1978. Dies sind (n-l)-mal stetig differenzierbare Funktionen, die auf den Gitterzellen 
mit Polynomen vom Grad n ubereinstimmen, wie in Fig, 4 dargestellt. Die obere Bildhalfte zeigt einen TVager Qk des B- 
Splines vom Grad n = 2, der Dimension d = 2 und der Glattheit m = 1 . In der unteren Bildhalfte ist der resultierende Ten- 30 
sorprodukt B-Spline bk abgebildet. Der TVager Qk, also die Vereinigung aller Gitterzellen, auf denen der B-Spline bk 
nicht identisch Null ist, besteht aus (n+l) d Gitterzellen; es gilt genauer Qk = kh + [0, (n+l)h] d . In alien Figuren ist der B- 
Spline bk an der Stelle kh markiert, also beispielsweise im Fall d = 2 in der linken unteren Ecke des TVagers. Fur die FE- 
S' imulation sind nur solche B-Splines von Belang, deren Trager das Simulationsgebiet Q schneiden; sie werden als rele- 
vante B-Splines bezeichnet. Die relevanten B-Splines werden nochmals in zwei Gruppen eingeteilt: diejenigen B-Spli- 35 
nes, bei denen der Teil des Tragers im Inneren des Simulationsgebietes groBcr als eine vorgegebene Schranke s ist, wer- 
den als innere B-Splines bezeichnet. Alle iibrigen relevanten B-Splines werden als auBere B-Splines bezeichnet. 
[0016] Die Aufgabe der vorliegenden Erfindung wird durch das im Anspruch 1 bestimmte Verfahren gelost. Besondere 
Ausfuhrungsarten der Erfindung sind in den Unteranspriichen bestimmt. Der Anspruch 11 bestimmt eine erfindungsge- 
maBe Computereinrichtung. 40 
[0017] In der rechten Bildhalfte der Fig. 1 ist die Einbindung des erfindungsgemaBen Verfahrens in den Ablauf einer 
FE-Simulation nach dem Stand der Technik bzw. die Substitution einiger Verfahrensschritte einer FE-Simulation nach 
dem Stand der Technik durch das erfindungsgemaBe Verfahren dargestellt. 

[0018] Die Eingabe 1 des Simulationsgebietes Q. kann iiber ein Eingabemittel erfolgen, insbesondere auch durch Ein- 
speichem von Daten, die aus einer computerunterstutzten Konstruktion (CAD/CAM) abgeleitet sind. Beispielsweise 45 
konnen so die bei der Konstruktion eines Kraftfahrzeugs verwendeten Daten unmittelbar in die erfindungsgemaBe FE-Si- 
mulation einflieBen. . 

[0019] Bei der Eingabe 2 und Abspeicherung des Typs der Randbedingungen wird zwischen natiirlichen und wesent- 
lichen Randbedingungen unterschieden. Die erfindungsgemaBe Basis wird fur homogene Randbedingungen gleichen 
T^ps konstruiert. Insbesondere verschwinden bei wesentlichen Randbedingungen die Basisfunktionen auf dem Rand T. 50 
Inhomogene Randbedingungen lassen sich dann beim Aufstellen des FE-Systems mittels Methoden, die dem Stand der 
Technik entsprechen, behandeln. 

[0020] SchlieBlich werden Steuerparameter eingclesen 3. Diese betreffen den Grad n und die Gitterweite h der zu vcr- 
wendenden B-Splines sowie die Schranke s zur Klassifikation der inneren und auBeren B-Splines. Bei fehlender Spezi- 
fikation konnen alle diese Eingabeparameter mit Hilfe von empirisch oder analytisch konstruierten Bewertungsfunktio- 55 
nen automatisch bestimmt werden. 

[0021] Die folgende erfindungsgemaBe Konstruktion der Basisfunktionen gliedert sich in die in Fig. 3 schematisch 
dargestellten Schritte, die nun beschrieben werden. 

[0022] Nach dem Einlesen des Simulationsgebietes Q. wird im ersten Verfahrensschritt ein das Simulationsgebiet £2 
uberdeckendes Gitter erstellt. Dann wird gepriift, welche der Gitterzellen ganz, teilweise oder nicht im Simulationsgebiet 60 
£1 liegen, die Zelltypen bestimmt 4 und diese Information iiber die Zelltypen gespeichert. Hierzu sind im Wesentlichen 
Innen/AuBen-Tests sowie Schnittbestimmungen zwischen dem Rand T des Simulationsgebietes Q und den Strecken 
bzw. Quadraten, die die Gitterzellen beranden, durchzufuhren. Fig. 7 zeigt die Ein- und Ausgabedaten fiir diesen Verfah- 
rensschritt. 

[0023] Im zweiten Verfahrensschritt werden zunachst anhand der Information iiber die Zelltypen die relevanten B- 65 
Splines bestimmt. Dann wird die Klassifikation 5 in innere und auBere B-Splines durchgefuhrt; die entsprechenden Li- 
sten der Indizes werden mit I und J bezeichnet. Hierzu ist mit Hilfe der im ersten \ferfahrensschritt gewonnenen Daten 
die GroBe derjenigen Teile der TVager der B-Splines zu beslimmen, die im Inneren des Simulationsgebietes liegen und 
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mit der vorgegebenen Schranke s zu vergleichen. Fig. 9 zeigt die Ein- und Ausgabedaten fur diesen Verfahrensschritt, 
[0024] Im dritten Verfahrensschritt werden Kopplungskoeffizienten ^ berechnet 6, die innere und auBere B-Splines 
gemaB der Vorschrift 

s Bi(x) = bi{x) + iei (i) 

verbinden. Jedem inneren B-Spline bi wird damit ein erweiterter B-Spline Bi zugeordnet. Die Konstruktion und die Ei- 
genschaften der Indexmengen J(i) und der Kopplungskoeffizienten eg sind wie folgt gegeben. Die Indexmengen J(i) be- 

10 stehen aus Indizes auBerer B-Splines. Ihnen entsprechen komplementare Indexmengen I(j) von Tndizes innerer B-Spli- 
nes; es gehort also i genau dann zu I(j), wenn j zu J(i) gehort. Fur einen gegebenen auBeren Index j ist die Indexmenge I(j) 
ein Feld, also ein quadratischer bzw. wiirfelfbrmiger Bereich von (n+l) d inneren Indizes, der durch einen minimalen Ab- 
stand zum Index j gekennzeichnet ist. FQr einen gegebenen auBeren Index j und einen inneren Index i in der Indexmenge 
I(j) sei pi dasjenige d-variate Polynom vom Grad n in jeder Variablen, das an der Stelle i den Wert 1 und an alien anderen 

15 Punkten des Feldes I(j) den Wert 0 hat. Dann ist der Kopplungskoeffizicnt als Wert von pi an der Stelle j gegeben, es 
gilt also ejj = p;(j). Die konkreten Werte der Kopplungskoeffizienten konnen entweder fur verschiedene Grade und rela- 
tive Positionen von j und I(j) tabelliert werden oder aber mit Hilfe von Lagrangepolynomen einfach berechnet werden. 
Fig. 1 1 zeigt die Ein- und Ausgabedaten fUr diesen Verfahrensschritt. 

[0025] Sind natiirliche Randbedingungen gegeben, so werden die in Gleichung (1) definierten erweiterten Splines ohne 
20 Modification zur weitcren Durchfiihrung des FE- Verfahrens verwendet. Sind dagegen wesentlichc Randbedingungen ge- 
geben, so ist noch eine Gewichtung gemaB der Vorschrift 

25 

durchzufuhren. Die zugehorige Abfrage erfolgt in einem option alen Verfahrensschritt 6a. Die so definierten Funktionen 
werden als gewichtete erweiterte B-Splines (WEB-Splines) bezeichnet. Formal gesehen entsprechen die bei naturlichen 
Randbedingungen verwendeten erweiterten B-Splines dem Sonderfall w(x) = 1. Sie werden deshalb ebenfalls als WEB- 
Splines bezeichnet. Fur den Fall wesentlicher Randbedingungen ist die Gewichtsfiinktion w wie folgt gekennzeichnet: 

30 FUr alle Punkte x des Simulationsgebietes lasst sich w(x) nach oben und unten durch von x unabhangige positive Kon- 
stanten mal den Abstand dist(x) des Punktes x vom Rand T abschatzen. Mit anderen Worten ist w im Inneren von CI po- 
sitiv und verschwindet in der Nahe des Randes T so schnell wie die Abstandsfunktion dist. Fiir Simulationsgebiete, die 
durch elementare geometrische Objekte (Kreise, Ebenen, Ellipsen, etc.) berandet sind, kann eine geeignete Gewichts- 
funktion gegebenenfalls in expliziter analytisc her Form angegeben werden. Anderenfalls sind Berechnungsvorschriften 

35 heranzuziehen, die typischerweise eine Glattung der Abstandsfunktion darstellen. Der Skalierungsfaktor l/w(xj entsteht 
durch Auswertung der Gewichtsfunktion im Gewichtspunkt xj. Dies kann ein beliebiger Punkt im TVager des B-Splines 
b, sein, der vom Rand mindestens den halben Schrankenwert s/2 entfernt ist. 

[0026] Als Resultat des erfindungsgemaBen Verfahrens erhalt man eine Berechnungs vorschrift fur die WEB-Splines Bi 
(vgl. Definitionen (1) und (2)), die alle giinstigen Eigenschaften gemaB der Aufgabenstellung aufweisen. Damit kann nun 
40 das FE-Verfahren gemaB dem Stand der Technik weitergefUhrt werden. Es ist dabei jedoch moglich, nach Anspruch 11 
die regulars Gitterstruktur der Basis funktionen vorteilhaft auszunutzen. 

[0027] Zusammenfassend ist zu sagen, dass fur das erfindungsgemaBe Verfahren unter anderem die Kopplung der au- 
Beren an die inneren B-Splines von Bedeutung ist. Diese bewirkt, dass die konstruierte Basis die fiir FE-Berechnungen 
wesentlichen Eigenschaften besitzt. Insbesondere ist eine erfindungsgemaBe Basis Bj (i aus der Indexmenge I) gleichma- 
45 Big in Bezug auf die Gitterweite h stabil und besitzt die gleiche Fehlerordnung wie die B-Splines b± bei Approximation 
glatter Funktionen, die die gleichen Randbedingungen erfullen. Zum anderen ist durch die Verwendung der Gewichts- 
funktion w die Einhaltung wesentlicher Randbedingungen gewahrleistet. 

[0028] Die zuvor beschriebenen Verfahrensschritte und deren Eigenschaften sowie die bevorzugte Ausfuhrungsart und 
die zugehorigen Ausfuhrungsbeispiele sollen mit Hilfe der nachfolgend aufgelisteten Figuren weiter erlautert bzw. illu- 
50 striert werden. Dabei konnen die in den Anspriichen und in der Beschreibung erwa^inten Merkmale jeweils einzeln fur 
sich oder in beliebiger Kombination erfindungswesentlich sein. 

[0029] Fig. 1 zeigt zum Stand der Technik die einzelnen Schritte beim Ablauf einer Finite-Elemente-Simulation und 
ordnet die erfindungsgemaBe Ermittlung der WEB-Basis in diesen Ablauf ein, 

[0030] Fig. 2a stellt einige Finite Elemente des Standes der Technik dem in Fig. 2b gezeigten WEB-Element gegen- 
55 iiber und listet die fur Finite-Elemente-Approximationen relevanten Parameter auf, 

[0031] Fig. 3 zeigt in einem Flussdiagramm die Verfahrensschritte zur Ermittlung der WEB-Basis, 
[0032] Fig. 4 zeigt einen TVager und den zugehorigen Tensorprodukt-B-Spline vom Grad 2, 

[0033] Fig. 5 illustriert die Problernstellung des ersten Ausfuhrungsbeispiels (Auslenkung einer Membran unter kon- 
stantem Druck) und zeigt die zugehorige Losung, 

60 [0034] Fig. 6 zeigt die Zelltypen anhand des ersten Ausfuhrungsbeispiels, 

[0035] Fig. 7 gibt eine ttbersicht der Ein- und Ausgabedaten des Verfahrens zur Bestimmung der Zelltypen, 
[0036] Fig. 8 illustriert exemplarisch die Klassifikation der B-Splines anhand des ersten Ausfuhrungsbeispiels, 
[0037] Fig. 9 gibt eine Obersicht der Ein- und Ausgabedaten des Verfahrens zur Klassifikation der B-Splines, 
[0038] Fig. 10 zeigt die Kopplungskoeffizienten eines auBeren B-Splines und die zugehorigen inneren B-Splines an- 

65 hand des ersten Ausfuhrungsbeispiels, 

[0039] Fig. 11 gibt eine Ubersicht der Ein- und Ausgabedaten des Verfahrens zur Berechnung der Kopplungskoeffi- 
zienten, 

[0040] Fig. 12 illustriert die Konstruktion der Gewichtsfunktion der bevorzugten Ausfuhrungsart, 
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10041] Fig, 13 zeigt den Trager eines WEB-Splines und die zugehorigen Kopplungskoeffizienten anhand des ersten 
Ausfiihrungsbeispiels, 

[0042] Fig. 14 erlautert die Problemstellung eines zweiten Ausfiihrungsbeispiels (inkompressible StrOmung) und de- 
ren Losung anhand der Stromungslinien und einer Abbildung zur Verteilung der Stromungsgeschwindigkeit, 
[0043] Fig. 15abis 15czcigenfurdaszwcitc Ausfuhrungsbeispicl zu den B-Spline-Gradcn n = 1,2,3 und gleicher Git- 5 
terweite die sich andernde Klassifikation der B-Splines, 

[0044] Fig. 1 6a bis 16c geben Auskunft Uber die Fehlerentwicklung bei der Finite-Elemente- Approximation rnit Hilfe 
von WEB-Splines sowie Uber das Laufzeitverhalten der WEB- Approximation fUr das zweite AusfUhrungsbeispiel, 
[0045] Fig. 17 stellt die WEB-Basis einem auf linearen Ansatzfunkdonen Uber einer TViangulierung (Stand der Tech- 
nik) basierenden Verfahren gegenuber und 10 
[0046] Fig. 18 zeigt eine erfindungsgemaBe Computereinrichtung. 

[0047] Eine besonders gunstige Ausfuhrungsart des erfindungsgemaBen Verfahrens, die als WEB-Verfahren bezeich- 
net wird, ist durch folgende Festlegungen gegeben: 

Die Schranke s wird so gewahlt, dass die inneren B-Splines dadurch gekennzeichnet sind, dass mindestens eine der Git- 
tcrzellcn ihres Tragers vollstandig im Simulation sgebiet ft liegt. Da zur Bestimmung der relcvanten B-Splincs ohnchin 15 
der Schnitt zwischen den Gitterzellen und dem Rand T zu berechnen ist, ist somit zur Klassifizierung kein wesentlicher 
zusatzlicher Rechenaufwand notwendig. Der Gewichtspunkt x^ wird als Mittelpunkt einer vollstandig im Simulauonsge- 
biet Q, liegenden Gitterzelle im Trager des B-Splines bi gewahlt. Auch dies ist effizient moglich, da die Bestimmung einer 
solchen Zelle bereits Bestandteil der Klassifizierungsroutine ist. 

[0048] Falls keine explizite analytische Form der Gewichtsfunktion bekannt ist, wird diese wie folgt angesetzt: 20 

1 falls dist(x) > 8 

w(x) (3) 

1 -(1 - dist(z)/J)" falls dist(x) < 6 25 

[0049] Fig. 12 illustriert die Konstruktion der Gewichtsfunktion. Dabei gibt der Parameter 5 die Breite eines Streifens 
Q$ an, innerhalb dessen die Gewichtsfunktion zwischen dem Wert 0 auf dem Rand des Simulationsgebietes und dem 
Wert 1 auf einem Plateau uber QJQ$ variiert. Der Parameter 5 wird so bestimmt, dass die Glattheit der Gewichtsfunktion 
gewahrleistet ist. 30 
[0050] Ein wesentlicher Vorteil des Verfahrens ist, dass keine Vemetzung des Simulationsgebietes erforderlich ist. 
Dies bewirkt in technischen Anwendungen eine deutliche Erspamis an Rechenzeit und Speicherkapazitat und verein- 
facht den Simulationsablauf. Dabei ist die Verfahrensstruktur fiir zwei- und dreidimensionale Probleme formal und auch 
in der technischen Durchfuhrung weitgehend identisch. Dies ermoglicht die zeit- und kostensparende Implementierung 
von Losern fur verschiedenste Anwendungen basierend auf einheitlichen Programmstrukturen. Die \ferwendung von B- 35 
Splines entspricht dem industriellen Standard bei der Modellierung geometrischer Objckte und schafFt damit eine natur- 
liche Verbindung zwischen FE- und CAD/CAM- Anwendungen. Bestehende umfangreiche Programmbibliothekcn aus 
beiden Bereichen konnen zur Implementierung einer auf dem erfindungsgemaBen \ferfahren beruhender FE-Simulation 
eingesetzt werden. Die nach dem WEB-Verfahren konstruierten Basisfunktionen besitzen alle Standardeigenschaften Fi- 
niter Elemente. Hierzu gehort insbesondere die Stabilitat der Basis. Diese bewirkt, dass beispielsweise bei linearen ellip- 40 
tischen Randwerlproblemen die Konditionszahl des entstehenden Gleichungssystems bei kleiner werdender Gitterweite 
nicht schneller anwachst als bei optimalen Triangulierungen. FUr die Anwendungen bedeutet dies beispielsweise, dass 
sich lineare Gleichungssysteme, wie sie bei der FE-Methode typischerweise entstehen, mittels iterauver \ferfahren effi- 
zient losen lassen. Weiterhin ist bei gegebenem Grad die Approximationsordnung maximal und die Anzahl der benotig- 
ten Parameter minimal. Damit sind sehr genaue Approximationen schon mit relativ geringer Parameterzahl moglich. 45 
Konkret kann dies bedeuten, dass Genauigkeiten, die bislang den Einsatz von GroBrechnem erforderten, auf Arbeits- 
platzrechnern erzielt werden konnen. Die regulare Gitterstruktur der erfindungsgemaBen Basis ermoglicht eine sehr effi- 
ziente Implementierung, insbesondere bei der Aufstellung und Losung von FE-Systemen. Ferner gelingt es durch \fer- 
wendung der Gewichtsfunktion Randbedingungen bei der Simulation einzuhalten, ohne die regulare Gitterstruktur der 
Basisfunktionen zu beeintrachtigen. SchlieBlich kann bei der Losung linearer elliptischer Randwertprobleme durch den 50 
Einsatz von Mehrgitter- Verfahren zur Losung des entstehenden linearen Gleichungssystems erreicht werden, dass die 
Gesamtlosungszeit proportional zur Anzahl der gesuchten Koeffizienten und damit optimal ist. 

[0051] Das erfindungsgemaBe Verfahren in der speziellen bevorzugten Ausfuhrungsart (WEB-Verfahren) wird anhand 
des in Fig. 5 gezeigten ersten Ausfiihrungsbeispiels veranschaulicht. Die Differenzialgleichung und Randbedingungen 
sind sehr elementar gewahlt, so dass ohne groBen Mehraufwand neben der erfindungsgemaBen Konstruktion der WEB- 55 
Basis der gesamte Ablauf der FE-Simulation verfolgt werden kann. 

[0052] In der oberen Bildhalfte ist eine elastische Membran gegeben, die entlang dem Rand T eines ebenen Simulati- 
onsgebietes £2 fixiert ist und auf die im Inneren des Gebietes ein konstanter Druck f = 1 wirkt. Die Verschiebung u erfullt 
bei geeigneter Normierung die Poisson-Gleichung mit homogenen Randbedingungen, 



60 



-Au = 1 in ft 
u =0aufT. 



[0053] Die Verschiebung u bzw. Auslenkung der Membran ist in der unteren Bildhalfte der Fig. 5 ersichtlich. Wie oben 
beschrieben, gliedert sich das WEB-Verfahren in die folgenden Schritte. 65 
[0054] Eingabe 1 des Simulationsgebietes Q: Der Rand T ist eine periodische Spline-Kurve vom Grad 6, die durch ihre 
Kontrollpunkte 20 (in Fig. 5 mit schwarzen Punkten gekennzeichnet) gespeichert wird. 

[0055] Eingabe 2 der Randbedingungen: Die homogene Randbedingung ist wesentlich, so dass die Konstruktion einer 
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Gewichtsfunktion erforderlich ist. 

[0056] Eingabe 3 der Steuerparameter: Es wird der Grad n = 2 und, urn die Figuren Ubersichtlich zu gestalten, eine re- 
lativ groBe Gitterweite h = 1/3 verwendet. 

[0057] Beslimmung 4 der Zelltypen: Wie in Fig. 6 illustriert, wird das Simulationsgebiet durch ein Gitter 21 Uberdeckt, 
5 das die Gittcrzcllen der TVagcr allcr potenzicll fur die Basiskonslruktion rclcvantcn B-Splincs cnthalt. Die Typbcstim- 
mung ergibt in dem Beispiel 69 auBere Gitterzellen 22 und 11 innere Gitterzellen 24, sowie 20 Gitterzellen 23 auf dern 
Rand. 

[0058] Klassifikation 5 der B-Splines: Der TVager des B-Splines bk ist hier das Quadrat Qk mit den Ecken 

10 (k b k 2 )h, (k x + 3, k 2 )h, (k, + 3, k 2 + 3)h, (ki, k 2 + 3)h 

[0059] Q(-4,o) und Q( 2 ,i) sind in Fig. 8 gezeigt. Die Gitterpunkte kh der relevanten B-Splines, fur die Qk das Innere des 
Simulationsgebietes schneidet, sind in Fig. 8 durch einen Punkt oder einen Kreis markiert. Alle Gitterpunkte ih flir innere 
B-Splines (i aus der Indexliste I), fur die mindestens eine Zelle des TYagers Qi vollstandig im Inneren von H liegt, sind 
15 durch cincn Punkt markiert. Beispiclswcisc liegt fiir i = (-4, 0) die Gittcrzelle (-2, 0)h + [0, h] 2 ganz in fit Alle Gitter- 
punkte jh fur auBere B-Splines (j aus der Indexliste J), fur die keine Zelle des TYagers Q ganz in Q. liegt, sind durch einen 
Kreis markiert. 

[0060] Berechnung 6 der Kopplungskoeffizienten: Zur Bestimmung der Kopplungskoeffizienten e^ wird fiir jedes fe- 
ste j der Indexliste J der nachstgelegene 3 x 3-Bereich 

20 

IC)=(li,li + Ui+2}x(l 2 ,l 2 +U 2 + 2} 

von Indizes in I gesucht. In Fig. 10 ist fur den mit einem Kreis markierten auBeren Gitterpunkt j = (-1, 2) der Bereich I(j) 
durch Punkte gekennzeichnet. Die Fig. 10 zeigt ebenfalls die zugehbrigen Kopplungskoeffizienten in Matrixdarstellung. 
25 Sie werden mit bivariater Interpolation berechnet. Bei spiels weise ist das Interpolationspolynom zu i = (-1, -1) 

p i (x) = -x l (x 1 +2)x 2 (x 2 -l)/2. 

[0061] Sein Wert an der Stelle x = j = (-1 , 2) ist pi(j) = 1 = ey, Es fallt auf, dass viele der Kopplungskoeffizienten 0 sind. 
30 Dies ist ein typisches Phanomen. Die Kopplungskoeffizienten eg sind nur dann fiir alle i der Indexliste I(j) ungleich 0, 
wenn die Indizes i in jeder Komponente von dem Index j verschieden sind. 

[0062] Berechnungsvorschrift 7 fur die Gewichtsfunktion: Die Gewichtsfunktion ist durch die Gleichung (3) gegeben 
mit n = 2 und 8 = 0.2. Der Parameter 6* wird numerisch berechnet. Er muss klein genug sein, so dass Singularitaten der 
Abstandsfunktion vermieden werden. Zur Berechnung der Abstandsfunktion erstellt das \ferfahren ein Programm, das 
35 mit Hilfe des Newton- Verfahrens arbeitet. Da die Gewichtsfunktion nur in einem Randstreifen ungleich 1 ist, ist der Auf- 
wand bei der spatcrcn Auswertung gering. 

Ausgabe 

40 [0063] Fig. 1 3 zeigt den TVager eines WEB-Splines B[ und die zu seiner Beschreibung benotigten Daten. Dies sind die 
Indexliste J(i) der mit b\ gekoppelten auBeren B-Splines bj, die Kopplungskoeffizienten e^ und der Gewichtspunkt X[. 
Diese Daten werden unter Einbeziehung der Gewichtsfunktion zur Erstcllung einer Berechnungsvorschrift fur die WEB- 
Splines benutzt. 

[0064] Der weitere Ablauf der FE-Simulation folgt dem Stand der Technik. 
45 [0065] Aufstellen 9 des FE-Systems: Die Eintrage der Systemmatrix und der rechten Seite sind 

G k ,i = Jngrad B k grad B if F k = J n f B k , k, i el. 

[0066] Das Gleichungssystem GC = F fiir die Basiskoeffizienten Ci hat im betrachteten Beispiel die Dimension 3 1 . Die 
50 Matrixeintrage G^i werden mit Hilfe von numerischer Integration berechnet, ebenso die Integrale F k . 

[0067] Losen 10 des FE-Systems: Das Galerkin- System wird iterativ mit dem Verfahren der konjugierten Gradienten 
gelost, wobei die SS OR- Vorkonditionierung zur Konvergenzbeschleunigung Verwendung findet. Nach 24 Iterations- 
schritten wird die Losung im Rahmen der Rechnergenauigkeit gefunden (Toleranz < le - 14). 

[0068] Berechnung 11 und Ausgabe 12 der Approximation: Die mit dem erfindungsgemaBen Verfahren berechnete 
55 Approximation ist u = liCiBi und wird in in der unteren Bildhalfte von Fig. 5 grafisch dargestellt. Der relative Fehler in 
der LrNorm betragt 0.028 . 

[0069] Die Leistungsfahigkeit des erfindungsgemaBen Verfahrens in der speziellen bevorzugten Ausfuhrungsart 
(WEB -Verfahren) wird in einem zweiten Ausfuhrungsbcispiel anhand der Simulation einer inkompressiblen Stromung 
veranschaulicht. Die in Fig. 14 gezeigte Anordnung zweier kreisformiger Hindernisse in einem parallel berandeten Ka- 
60 nal dient der Illustration des prinzipiellen Vorgehens. Bei komplizierteren Geometrien, wie sie in konkreten Anwendun- 
gen typischerweise vorliegen, arbeitet das Verfahren volikommen analog und effizient. In Fig. 14 sind die Stromungsli- 
nien 25 innerhalb des durch Ti bis T 4 sowie durch T 5 und T 6 berandeten Gebiets dargestellt. Die Differenzialgleichung ist 

Au = 0 in CI 

65 

mit den Randbedingungen 
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du du du 

— = v Q auf I\ , — = -v 0 auf T 2 , ^ = 0 auf r 3) . . . , T 6 . 

[0070] Die Flussgeschwindigkeit v = -gradu ist in dcr untercn Bildhalfte daigcstellL 5 
[0071] Die erfindungsgemaBe Konstruklion der WEB-Basis lauft vollkommen analog zu dem ersten Ausftthrungsbei- 
spiel ab. Der einzige Unterschied ist, dass aufgrund der natiirlichen Randbedingungen keine Gewichtsfunktion erforder- 
lich ist. 

[0072] Fig. 15a bis Fig. 15c zeigen die Klassifikation der relevanten B-Splines fiir verschiedene Grade n (s. auch Fig. 
8). In der Abbildung sind die inneren B-Splines bj, die in die WEB-Basis ohne Erweilerung iibemommen werden, durch 10 
schwarze Dreiecke markiert. Fiir kleines h nirnmt die Anzahl dieser B-Splines zu, d. h., es gilt B\ = b, fur den uberwie- 
genden Teil der WEB-Basis. Im betrachteten Beispiel trifft dies bei Grad n = 3 f ur 236 von 252 Indizes i zu. 
[0073] Fig. 16a zeigt in zwei Diagrammen den numerisch bestimmten relativen l^-Fehler des Potenzials (linke Bild- 
halfte) in Abhangigkeit von der Gitterweite h = 2" k mit k = 1, . . ., 5 und die numerisch geschatzte Konvergenzordnung m 
(rechte Bildhalfte). Dabei werden fiir die verschicdenen Grade dcr WEB-Splines als Kennzeichncr * (n = 1), O (n = 2), 15 
A (n = 3), □ (n = 4) und I (n = 5) verwendet. ErwartungsgemaB ergibt sich m « n + 1, d. h. eine ungefahre Fehlerre- 
duktion um den Faktor 2 n+l bei Halbierung der Gitterweite. Analog ergibt sich fur den in der Fig. 1 6b abgebildeten rela- 
tiven Approximationsfehler der Flussgeschwindigkeit (H l -Norm der L6sung, linke Bildhalfte) eine Konvergenzordnung 
von m » n (rechte Bildhalfte) und die damit verbundene Fehlerredukuon um etwa den Faktor 2 n bei Halbierung der Git- 
terweite. 20 
[0074] Fig. 16c (rechte Bildhalfte) zeigt die Rechcnzeit in Sekunden fur die Konstruktion der WEB-Basis in Abhan- 
gigkeit von der Zahl der resultierenden Basisfunktionen gemessen auf einem Pentium H-Prozessor mit 400 MHz. Bei- 
spielsweise werden fiir die Konstruktion einer WEB-Basis vom Grad 3 zur Gitterweite h = 0.125 mit 2726 WEB-Splines 
1 .32 Sekunden benotigt. Es fallt auf, dass der Aufwand zum Erstellen der WEB-Basis weitgehend unabhangig vom Grad 
n der Basis ist. In der linken Bildhalfte der Fig. 1 6c ist die Anzahl der CG-Iterationen bezogen auf die Anzahl der Basis- 25 
funktionen gezeigt. So werden fur das zugehorige System mit 2726 Unbekannten 65 PCG-Iterationcn benotigt. Die Gc- 
samtrechenzeit inklusive Aufstellen und Losen des Galerkin-System betragt etwa 2.48 Sekunden. 
[0075] In Fig. 17 wird die WEB-Methode mit einem Standardlosungsverfahren, das das Simulationsgebiet vernetzt 
bzw. trianguliert (obere Bildhalfte) und Hutfunktionen verwendet, verglichen. Die Grafik zeigt in der unteren Bildhalfte 
den L 2 -Fehler bezogen auf die Anzahl der Parameter, Die Ergebnisse des Standardlbsers sind mit fetten Karo-Kennzeich- 30 
nern markiert und werden den mit Hilfe von WEB-Basen der Grade 1 bis 5 erzielten Ergebnissen gegenubergestellt. Bei- 
spielsweise wird eine Genauigkeit von 10~ 2 mit der WEB-Methode durch Verwendung von 213 Basisfunktionen mit 
Grad 2 erreicht bei einer Gesamtrechenzeit von 0.6 Sekunden. Zum Erzielen der gleichen Genauigkeit benotigt das Stan- 
dardverfahren mit linearen Hutfunktionen 6657 Basisfunktionen. 

[0076] Bei der Beurteilung des Standardlosungsverfahrens sind noch zwei weitere Aspekte zu beachten. Zum einen 35 
verdeutlicht Fig. 17, dass selbst eine moderate Genauigkeit von 10~ 3 mittels Hutfunktionen nur erreichbar ist, wenn weit 
mehr als eine Million Koeffizienten verwendet werden. Dies zeigt, dass genaue Resultate bei der Verwendung von Hut- 
funktionen in der Regel eine enorme Rechen- und Speicherkapazitat erfordern oder aber nach dem Stand der Technik gar 
nicht erreichbar sind. Zum anderen steigt der zur Vemetzung erforderliche Aufwand mit der Komplexitat des Simulati- 
onsgebietes. Im Gegensatz zu realistischen Anwendungen ist das hier untersuchte Gebiet aufgrund seiner sehr einfachen 40 
Struktur noch vergleichsweise einfach zu triangulieren. 

[0077] Das zweidimensionale Beispiel zeigt den Performancegewinn durch das WEB-Verfahren. Eine noch groBere 
Leistungssteigerung ist bei dreidimensionalen Problemen moglich. Zum einen ist dort der Aufwand fur die \fernetzung, 
die bei dem WEB-Verfahren entfallt, wesentlich groBer. Zum anderen macht sich hier auch die Reduktion der Anzahl der 
benotigten Basisfunktionen noch wesentlich deutlicher bemerkbar als im zweidimensionalen Fall, 45 
[0078] Die Fig. 18 zeigt eine erfindungsgemaBe Vorrichtung, insbesondere eine Computereinrichtung 30, mit Einga- 
bemitteln 31, 32, 33, Ausgabemitteln 34 sowie eine den Verfahrensablauf steuemde Steuereinrichtung 35. Zur Abarbei- 
tung des erfindungsgemaBen Verfahrens und insbesondere zum Zwecke der Parallelisierung des zugehorigen Rechenab- 
laufs bedient sich die zentrale Steuereinrichtung 35 vorzugsweise mehrerer arithmetischer logischer Einheiten (Arithme- 
tic Logical Unit, ALU) oder sogar mehrerer zentraler Verarbeitungseinheiten (Central Processing Unit, CPU) 36. Diese 50 
erlauben insbesondere eine parallele Abarbeitung der Verfahrensschritte Klassifizieren 5 der B-Splines, insbesondere 
auch Schneiden des regularen Gitters mit dem Simulationsgebiet ft, Ermitteln 6 der Kopplungskoeffizienten eg und/oder 
Auswerten der Gewichtsfunktion w(x) an Punkten x des Simulationsgebietes ft. 

[0079] Die Recheneinheiten 36 greifen dabei auf den gemeinsamen Datenbestand von Speichermitteln 37 zu. Die Ein- 
gabe der Daten kann beispielsweise mittels einer Tastatur 31, einem maschinenlesbaren Datentrager 38 iiber eine zuge- 55 
horige Lesestation 32 und/oder iiber ein drahtgebundenes oder drahtloses Datennetz mit einer Empfangerstation 33 er- 
folgen. iiber die Lesestation 32 bzw. einen zugehorigen Datentrager 38 kann auch das den Verfahrensablauf steuernde 
Steuerprogramm eingegeben werden und beispielsweise auf dem Speichermittel 37 dauerhaft abgelegt werden. Entspre- 
chend kann das Ausgabemittel 34 einen Drucker, einen Monitor, eine Schreibstation fiir einen maschinenlesbaren Daten- 
trager und/oder eine Sendestation eines drahtgebundenen oder drahtlosen Datennetzes aufweisen. 60 

PatentansprUche 

1. Verfahren zur Erhohung der Leistungsfahigkeit einer Computereinrichtung zur Finiten Elemente-Simulauon 
durch automatische Generierung geeigneter Basisfunktionen mit Hilfe von uniformen Tensorprodukt B-Splines mit 65 
den folgenden Schritten: 

- Definieren (1) eines Simulationsgebietes (12) und Abspeichern von Gebietsdaten des Simulationsgebietes 
(ft), 
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- Eingabe (2) und Abspeichern von Randbedingungen, 

- Fesllegen (3) einer vorgebbaren Gitterweite (h) und eines vorgebbaren Grades (n) der B-Splines, 

- Ermitteln (4) eines sich iiber das Simulationsgebiet (ft) erstreckenden Gitters und eines Typs der Gitterzel- 
len, 

5 - Klassifizicren (5) der B-Splincs, dcren Tragcr das Simulationsgebiet (ft) schneiden, durch Ermitteln von in- 

neren und auBeren B-Splines, wobei fur innere B-Splines der Schnitt des Tragers mit dem Simulationsgebiet 
(ft) ein Quadrat bzw. einen WUrfel einer vorgebbaren Kantenlange (s), welche mit der Weite des Gitters (h) 
ubereinstimmen kann, enthait, 

- Ermitteln (6) von Kopplungskoeffizienten (ey) zur Bildung von Linearkombinationen von inneren und au- 
10 Beren B-Splines durch Lagrange-Interpol ation auf zu den auBeren Indizes nachstliegenden Feldern innerer In- 

dizes, 

- Speicherung und Ausgabe der die Basisfunktionen bestimmenden Parameter. 

2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass vor dem Speichem und Ausgeben der Parameter der 
Schritt durchgefuhrt wird: 

15 Fcstlegen (7) einer vorgebbaren Gewichtsfunktion (w) zur Einhaltung wescntlicher Randbedingungen sowic Be- 

stimmen (8) von Gewichtspunkten und Skalierungsfaktoren. 

3. Verfahren nach Anspruch 2, dadurch gekennzeichnet, dass die Gewichtsfunktion w durch einen glatten Ober- 
gang von einem konstanten Plateau im Inneren des Simulationsgebietes (ft) zum Wert Null auf dem Rand (0 fest- 
gelegt wird. 

20 4. Verfahren nach einem der Anspruche 1 bis 3, dadurch gekennzeichnet, dass B-Splincs mit mindestens einer ganz 

im Simulationsgebiet (ft) enthaltenen Gitterzelle des Tragers als innere B-Splines klassifiziert wcrden. 
5. Verfahren nach einem der Anspruche 1 bis 4, dadurch gekennzeichnet, dass der Gewichtspunkt als Mittelpunkt 
einer ganz im Simulationsgebiet (ft) enthaltenen Gitterzelle des Tragers des entsprechenden inneren B-Splines ge- 
wahlt wird. 

25 6. Verfahren nach einem der Anspruche 1 bis 5, dadurch gekennzeichnet, dass das Simulationsgebiet (ft) durch 

Einspeichern von aus einer computergestutztcn Konstruktion (CAD/CAM) ableitbarcn Daten definiert wird. 
7. Verfahren nach einem der Anspruche 1 bis 6, dadurch gekennzeichnet, dass die Gitterweite h unter \ferwendung 
von abgespeicherten empirisch und/oder analytisch gewonnenen Werten mittels einer zugehorigen ersten Bewer- 
tungsfunktion automatisch festgelegt wird. 

30 8. Verfahren nach einem der Anspruche 1 bis 7, dadurch gekennzeichnet, dass der Grad n unter Verwendung von 

abgespeicherten empirisch und/oder analytisch gewonnenen Werten mittels einer zugehorigen zweiten Bewertungs- 
funktion automatisch festgelegt wird. 

9. Verfahren nach einem der Anspruche 1 bis 8, gekennzeichnet durch die Schritte: 

- Aufstellen (9) eines bei der FE-Simulation zu losenden Gleichungssystems, 
35 - Losen (10) des Gleichungssystems, 

- Berechnen (11) einer Naherungslosung, und 

- Ausgeben (12) der Naherungslosung. 

10. Verfahren nach Anspruch 9, dadurch gekennzeichnet, dass beim Losen (10) des Gleichungssystems ein Mehr- 
gitter- Verfahren eingesetzt wird. 

40 11. Vorrichtung zur Durchfuhrung eines Verfahrens nach einem der Anspruche 1 bis 10, insbesondere Computer- 

einrichtung, mit Eingabemitteln (31, 32, 33) und Ausgabemitteln (34), Speichermitteln (37) und einer Rechenein- 
heit (35, 36), dadurch gekennzeichnet, dass die regulare Gittcrstruktur zur Optimierung des Rechenablaufs, insbe- 
sondere durch Parallelisierung, ausgenutzt wird. 

12. Maschinenlesbarer Datentrager (38), insbesondere Magnetband, Magnetscheibe, Compact Disc (CD) oder Di- 
45 gital Versatile Disc (DVD), dadurch gekennzeichnet, dass auf dem Datentrager (38) ein Steuerprogramm fur eine 

Computereinrichtung (30) gespeichert ist, gemaB dem durch die Computereinrichtung (30) ein Verfahren nach ei- 
nem der Anspruche 1 bis 10 abgearbeitet werden kann. 
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Abstract of DE1 0023377 

A simulation area ( OMICRON ) is defined and 
the data are stored and marginal conditions are 
input and stored. Grid area and grid cells 
extending over the simulation area are 
determined, to determine and classify inner and 
outer B-splines of preset grid width and degree 
(n) in the simulation area. Coupling coefficients 
(eij) for linear combinations of inner and outer B- 
splines are evaluated, which are stored and 
output as parameters for the basic functions. An 
Independent claim is also included for computer 
device with improved finite element simulation 
based basic functions. 
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